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TITL OTSerQuRs - REAL ** INTEGER*4 power routine 
»IDENT /1-006/ ; File: OTSPOWRJ.MAR Edit: $BL1006 


’ 
FRRAAKAK SHH KAH HEHE T AHHH HTH HAHRAHKAARRAHHHAKeHKeeHKeKeeeeeeeteeeeeeeeeeeeeee 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY He ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION oe THE ao gid NOTICE. THIS SOFTWARE OR ANY OTHER 


Vv 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


;* 

‘@ a 
*¢ * 
*@ ® 
*@ o 
*@ a 
*@ e 
*®@ 
*@ * 
*@ 
;* TRANSFERRED. : 
** * 
*@ * 
*® ® 
*@ ® 
*@ = 
*@ ® 
*@ ® 
*® * 
*@¢ * 
*@ a 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
CoRPORAT Toke NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


ARBAB RRRRRRARARRRALAASESAS SEER SAREE SESE LER ER EERE EERE EERE SEE SESE R ECE E SEES 


s FACILITY: Language support Library - user callable 
; ABSTRACT: 


Floating base to integer senguere power. 
Floating overflow and underflow can occur. 


HISTORY: 
AUTHOR: 
Thomas N. Hastings, 5-May-77: Version 0 


; modified by: SUSAN HUBBARD AZIBERT 


ODNA NEW OB OD NA NE WIN OOD NOUN EWI OC OO NOUN EW "OOOO 


: VERSION: 0 


o V04-00 Page 


Undefined exponentation can occur if base is 0 and power is 0 or negative. 


fF ? | 

REAL ** INTEGER®4 power routine i ob dat BS 9 96:03:88 AX/VMS Macro V04-00 Page § 
HISTORY ; Detailed Current Edit. History 6=SEP=1984 11:28: MTHRTL. SRCJOTSPOWRJ MAR; 1 (3) | 
1 .SBTTL HISTORY ; Detailed Current Edit History 


version = changed module name to forpow 
version = changed error handler from MTNSERROR to MTHSSERROR 
version 08 = removed W* from MTHSSERROR, saved code with MOVZBL. 
removed infinite loop with largest negative laneae exponent. 


e § on 09 - changed MTHSSERROR to MTHSSSIGNAL = JMT 
- - Update version number and copyri ht no ice. The previous 
version number was 0-10. J 
change, Pao UNDEXP to MTHSK UNDEXP. , 85, 507 -DEC- -78 
Add o the PSECT of H5tive:., 
Declare Rind yy mn? 
Add handlers fe Steck §s$ FLTOVF and SS$_FLTDIV., and signal 
MTHS_FLOOVEMAT or MTHS_FLOUNDMAT instead, depending on the context. 
JAW 26-Feb-1980. 
1-006 = Use general mode addressing. SBL 30-Nov-1981 


Edit ae for Version 0 of OTSSPOWRJ 
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-SBTTL DECLARATIONS 


INCLUDE FILES: 


; EXTERNAL SYMBOLS: 


DSABL 

sEXTRN MTHSK_UNDEXP, MTHSK -FLOOVEMAT MTHSK har Ss 

~EXTRN BTHSSSIGNAL Math error routin 

EXTRN SSS$_FLTOVF, SS$_FLTOVF_F, ses -FLTDIV, SS$_ PLIDIV. F, SS$_CONTINUE 
> MACROS 
‘ SCHFDEF ; Define condition handler symbols. 

SSFDEF ; Define stack frame symbols. 

SPSLDEF ; Define program status longword 


symbols. 


: EQUATED SYMBOLS: 


base = 4 
exp = 


base input formal - by-value 
exponent intpu formal - by-value 


> OWN STORAGE: 


: PSECT DECLARATIONS: 


-PSECT _OTSSCODE PIC,SHR,LONG, EXE, NOWRT 
; program section for OTS$ code 


H 7 
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-SBTTL OTSSPOWRJ = floating to power longword giving floating result 


+e 
FUNCTIONAL DESCRIPTION: 


Floating result = floating base ** signed longword exponent 
The floating result is given by: 


base exponent result 
any >0 product (base * 2**i) where i is each 
non-zero bit position in exponent 


= 0 1.0 
. 8 ergetraes exponentation 


v Atv 
o ooo 


<0 1.0 / product (base * 2**i) 

where i is each non-zero bit position 
in ‘exponent: 

Undefined exponentation 

1.0 / product (base * 2**i) 


Al 

oo 
AA 
oo 


where i is each non-zero bit position 
in texponent: 


; Floating overflow can occur on either of the two MULF's. If this 

3 papeene when the exponent is less than zero, the ongept cen is caught by 
3 a local condition handler named EXC_HNDLR_UNDER, which sets the result 
3 to 0.0 and either signals MTH$_FLOUNDMAT Tif FU is enabled in the 

3 caller's PSW) or continues at POWRJX. If it happens when the exponent 
3 ‘s groomer than zero, the exception is caught by a local condition 
3 an 


COOOCCCCOCOCOCOCOOOCOOOOCOOOOOOOOOCoOO 


COOoOooOoCoooooooooooooO 
SOoOoOoOoCoCoOOOoOOoOOoOoOoOoOOOOoOoSoO 
SOOOoOoOOooOoooooooooooo 


ler named EXC_HNDLR_OVER, which sets the result to the reserved 
operand (-0.0) and signals MTH$_FLOOVEMAT. 


Floating overflow and floating divide by zero can occur on the DIVF. 
These exceptions are caught by EXC_HNDLR_ OVER, which sets the result to 
the reserved operand (-0.0) and signals ATHS_FLOOVEMAT. 


Undefined exponentiation occurs if base is 0 and 
exponent is 0 or negative. 


CALLING SEQUENCE: 
Power.wf.v = OTSSPOWRJ (base.rf.v, exponent.rl.v) 
INPUT PARAMETERS: 
NONE 
; IMPLICIT INPUTS: 
The setting of FU in the caller's PSwW. 
OUTPUT PARAMETERS: 
NONE 


Oo 

So 

o 

oO 
a cee ee ed od od od 
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DNA NEW OS OD NAME WIN O OD NA UNE WIN CO OD NOUN EWN O OD NOU EWI OO OONOUE WP 
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; IMPLICIT OUTPUTS: 
NONE 


~~ 


Sd 


Me 
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FUNCTION VALUE: | 
| 

Floating base ** exponent 
SIDE EFFECTS: 


Signals MTHS FLOOVEMAT if f ae Fg) pbs occurs on either of the two 
e S when X. onent > f floating overflow or divide by zero 


IVF. 
Signals s nine FL FLOUNDAAT ty floating overflow occurs on either of the two 
ULF's when gr onent and cailer has FU enabled. 
Signals MTHS UNDEXP (82 = “, OURD EE NED exponentation’) if 
base is 0 — exponent is 0 or negative. 


0004 ENTRY OTSSPOWRJ, “M<R2> disable integer overflow 
gecure on largest negative integer exp 
6D 66°AF 9E MOVAB B*EXC_HNDLR_OVER, (FP) anslate exceptions to 
MTHS F FLOOVEMAT. 
50 50 MOVF #1, RO RO = “initial result 
51 04 aC 50 MOVF base(AP), R1 R1 = base 
52 O08 AC 00 MOVL exp(AP), R2 R2 = exponent 


Note: integer overflow can occur 

on largest neagtive integer exponent 
However, R2 is correct unsigned 32-bit val 
Use ROTL -1 rather than ASHL -1 below. 
branch if exponent > 0 

Translate exceptions to 

MTH$_FLOUNDMAT. 


BGTR EX 


1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 PGTR 

1 MOVAB B“EXC_HNDLR_UNDER, (FP) 


0D 
6D S7*AF GE 


; test b 


51 33 R1 3 ase 
¢§ 1 BEQL UNDEF INED ; undefined 0**0 or 0**#(-n) 
52 2 cf MNEGL R2, R2 ; R2 = texponent: 
26Ci«d1 EQL POWRJX : if exponent is 0, return RO = 1.0 


s¢ 
: Exponent is > 0 or (exponent is =< 0 and base is not = 0 -=- use ‘exponent:) 


EXPGTR: BBSC #0, R2, PARTIAL 


SQUAR: ROTL #-1, R2, R2 
SQUAR1: MULF R1, R1 


08 52 00 €4 
m wow HH 9¢ 
51 1 


branch if itexponent: is odd 
ane clear low order bit 
= “bit unsigned exponent: /2 
Rf = current power iH Pee 
Floating atthe wi ree r. Toute 
and signal SS$_FLTO VE ‘3 SS$_F 
branch if next™bit in Mn ‘ -° 


SOOO SSS SSS SSS SSSSsssssss 
oo OOCCCOCOCCOOSSSOSOSSSSCSCOSOOOOOSOSOSOOOOOOOOOOSOOOSOoO 


POR OROPOPORIPONINIPITONIIITITIPININIIND — — — 9 3 3 3 3 33 3 SQOODOOCOOCOOCOOOCOOOCOSOOOSOOSSOO 


TRH HAHAHA RAIMAOLLOCCCOCOMD ON NNW 2 SF OPO ANNO OOCOCOCOOOCOCSOO 


BLBC R2, SQUAR 


7 Here when bit i of ‘exponent: is a 1 
: Partial result = partial result * ( 


PARTIAL: ; 
MULF R1, RO ; RO = new partial result 


7 
7 
7 
7 
7 
7 
7 
7 
2 
9 
g 
9 
G 
9 
y 
9 
9 
0 
0 
0 
0 
0 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
: base * 2**j) 


MEW @ OC ODA UNE WP 9 OD NA NE WWIN 9 OD NAME WIN 9 OD NAME WN 0 ODNAUES W000 


MTHS$_FLOOVEMAT. 
RO = 1.0/result 
return, result in RO 


50. 0608 «= 550s 47 DIVFS RO, #1, RO 


POWRJX: RET 


J? | 
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$ ASHL #-1, R2, R2 3; R2 = exponent! /2 | 

BNEG  SQuUARI ; loopback if more exponent bits are 1 

TSTL exp (AP) est sign of exponent | 

BGTR POWRJX f exponent > 0, return RO 

MOVAB B“EXC_HNDLR_OVER, (FP) Translate exceptions to 


g¢ 
3; Undefined exponentation error - 0**0 or 0**(-n) 


PPE LEE EE ES 


UNDEF INED: 
ASHL #15, #1, RO 

MOVZBL M#MTHSK_UNDEXP, -(SP) 

CALLS #1, G*ATHSSSIGNAL 


RO = reserved floating operand 
Indicate unde t toed exponentiation. 
convert to 32-bit condition code 


or 78 ; 
F OA ; 
; and SIGNAL MTH$_UNDEXP 


52 S52 FF BF 7 
aim 
08 a dS 
0 14 
6D 66°AF  9E 
0 
0 
50 01 0 
43 00° 0 
00000000'GF 01 FB 
; Note: 2nd arg not needed since 
no JSB OTS$POWRJ 


RET return 


+ 
The following handler is established to process exceptions which imply 
underflow of the final result (floating overflow in either of the two MULF's 
when exp < 0). On the occurrence of such an exception, the handler signals 
MTHS$_FLOUNDMAT. 


EXC_HNDLR_UNDER: 
-OORD “M<R2, R3, R4> ; Entry mask ; 
BSBB SETUP ; Set up RO:R3 and identify condition. 
;_ Return only if FLTOVF or FLTDIV. 
BBC #PSLSV_FU, SFSW_SAVE_PSW(R2), CON_U 
: Branch Tf caller has not enabled FU. 
even. eee RAAT R4 ; Report MTHS$_FLOUNDMAT, not SSS$_FLTOVF. 


28 
1€ 04 A2 06 1 
54 O0°SF 9A 
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PAEAAAAAA AA MAMMA MMA TUT 
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POOAOAA Ff OO BD ONIN OA A 9D NNN INIOR MO OOTY 


ROOD NA NE WR @OOONA UE WIN CO OO NAW EWN 0 OD NOU EWN ( O OONOU ES Wi OOo 


000909 NINN NIN NEN NINO AAA AAA AONUMA IIE BE BRR RRS 


0c 11 
s@ 
0 3; The fol lowing handler is established to process exceptions which en gy 
0 3; overflow of the final result (floating overflow in either of the two MULF's 
0 3 when exp > 0, floating overflow in the DIVF, or speet ine divide by zero in the 
bee > DIVE). On the occurrence of such an exception, the handler signals 
; MTHS_FLOOVEMAT. 
bee 7° 
6 EXC_HNDLR_OVER: 
0or¢ re -@ORD “M<R2, R3, R4> 3 Entry mask 
—.. 3 6 BSBB SETUP ; Set up RO:R3 and identify condition. 
6A ; Return only if FLTOVF or FLIDIV. 
50 01 of 78 6A ASHL #15, #1, RO ; Make the default result -0.0. 
54 00°8F 9A $ rOVZBL #MTHSK_FLOOVEMAT, R4 ; Report MTHS_FLOOVEMAT, not SSS_FLTxxx. 
10 A2. DD 7? DO_SIG: PUSHL SFSL_SAVE_PC(R2) ; Report caller's PC, not exception PC. 
33 DD 7? PUSHL R4 ; Report MTH$_xxx, not SS$_xxx. 
00000000'GF 02 FB 007 CALLS #2, G*MTHSSSIGNAL : Signal the Condition. 


K 7 
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CON_U: MOVQ RO, CHFSL MCH SAVR (R3) ; If continued, restore RO and R 
OV S*#SS$_ Con INDE. R ; Continue from the original F aettien 


M 
DO_RET: RET i; Exit from handler. 


In the case of a normal A. OF, (FLTOVE, FLTOVF _F, 


i equivalent to ceitut RES 
‘ t ge sets up RO:R3 as follows: 


FLTDIV 


TEG 
flo 
7 3 
4 
fis 
re 
8 ; Common setu cpgt ne for penglers, Returns normally if exception was FLTOVF, 
: Fonovr uF. FLIDI ert TOIV If the exception uae a nngta'e ng else, it 
30 > executes a PRET ing an an ai it from the handler with which is 
: 35 3 or FLTBIV. F) 
93 ; RO/R1: 
94; Re: Lstoons of establisher's frame 
é $6 3 R35: address of mechanism array 
6 297° 
50 7C Hf 98 SETUP: CLRQ RO ; ist default result to 0.0. 
52 O04 AC 7D 4 99 MOVa CHFSL_SIGARGLST(AP), R2 ; = address of signal array 
C 00 = , seerese of mechanism array 
0000'8F 04 a2 81 008 01 CMPW = CHF S$L_SIG_NAME(R2), wsss’ -FLTOVE 
9 8 ; as it a floating overflow trap? 
. ae 9 0 BEQL DO_RSB 3 pe if yes. 
0000 * 8F 04 A261 a 04 CMP CHFSL_SIG_NAME(R2), #SSS$_FLTOVF_F 
A 05 ; Ora Floating overflow fault? 
10 13 9A BS BEQL DO_RSB Branch if yes. 
0000'8F 04 A2 B81 009 0 CMPW  CHFSL_SIG_NAME(R2), #SS$_FLTDIV 
OA 08 : Or a floating divide by zero trap? 
08 13 OOA 09 BEQL DO_RSB Branch if yes. 
0000 ' 8F 04 A2 Bi OO0A4 10 CMPW CHFSL _SIG_NAME(R2), #SSS. _FLTDIV_F 
ae 11 ; Or a Float ting divide by zero fault? 
D9 Ss «12 ne \¢ BNEQ DO_RET : None of the above: return from handler 
08 A2 97 AF 9E boat 14 DO_RSB: MOVAB B*POWRJX, CHFSL_SIG NAMES 4 (RE) 
0081 15 ; Change return PC to POWRJX. 
52 04 a3 00 aed 16 MOVL CHFSL_MCH_FRAME(R3), R2 ; R2 = address of establisher's frame 
05 0085 17 RSB : Return. 
0086 18 
00B6 319 - END 


L 
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OTSSPOWRJ 
Symbol table 


1G 
XC_HNDLR_OVER 
XC"HNDLR_ -UNDER _ 


GTR 
MTHSSSIGNAL etereree x 
MTHSK_FLOOVEMAT ‘teeeeeee x 
x 


MTHSK_FLOUNDMAT teeeeeee 
hy 9 UNDE XP eeeeeene 


S$POWRJ eto 44 RG 
PARTIAL OOOOO2F R 

RJX 0000046 R 
PSL$V_FU = 00000006 

TUP iets 3 R 02 
SFSL_SAVE_P = 0000001 
SFSW SAVE _PSW = 444 64 4 

UAR 0000024 R 02 
SQUAR1 00000029 R Be 
SS$_CONTINUE eeeeeeee 8 6X 
SS$_FLTDIV seceneee§ =X 00 
SSS_FLIDIV_F eteceree =X 600 
SS$_FLTOVF eeeeecee xX 00 
SS$_FLTOVE _F eeceeeee §6X = 600 
UNDEF INED 00000047 R 02 

des on eer ereee eee eae + 
: ; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
. ABS . 00000000 < 9} NOPIC USR CON S LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS 00000000 < Q:) -) NOPIC USR CON ABS LCL NOSHR EXE RD — WRT NOVEC BYTE 
-OTSSCODE 00000086 ( 182.) 2°} PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
ey 
' Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
00.96 


Initialization 30 :00: 9 
Command processing 1 :00: 
ee tabl g 5 : ; 
adie sor :00:Q00. 
Pass 2 7% :00:00. 
Symbol table output $ :00:00. 
0 0:00: 


Psect synopsis output 
Cross-reference output 


04-00 Page 
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VAX=-11 Macro Run Statistics 
Assembler run totals 361 00:00:03.47 00:00:14.67 


ibs working set Limit was 900 pages 
8 bytes (17 pages) of virtual memory were used to buffer the intermediate code. 
* re vere 10 pages of symbol table ,feaee al ty d to hold 106 non-local ane local symbols. 


319 source Lines were read in Pass 1, produc object records in Pass 
pages of virtual memory were = to defi ne 9, macros. 


tow ate wre ee eee ee ee eee =-+ 


! Macro Library statistics ! 


werent ere ene ee + hth DA te | 


Macro Library name Macros defined 
“S255SDUAZB:(SYSLIBISTARLET.MLB;2 0 6 
148 GETS were required to define 6 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWRJ/OBJ=OBJ$:OTSPOWRJ MSRC$:O0TSPOWRJ/UPDATE=(ENHS$:OTSPOWRJ) 


T 
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